Dispersion anomalies induced by the low-energy plasmon in the cuprates 
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We discuss the characteristic effects of the eiectron piasmon interaction resulting from the ~ 
1 eV plasmon, which is a universal feature in the cuprates. Using the framework of a one-band 
tight binding model, we identify signatures of this low energy plasmon in the electronic structure of 
metallic overdoped Bi2212 as well as half-filled insulating SCOC. The electron-plasmon interaction 
is found to yield renormalizations near the Fermi energy in reasonable accord with experimental 
observations, and to produce dispersion anomalies at higher energies. 
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Quasiparticle dispersions near the Fermi energy in the 
cuprates are renormalized (reduced) by a factor of Z ~ 
0.3 — 0.6 [1] in comparison to the band theory predictions 
based on the conventional LDA picture, interrupted by 
the appearance of low energy 'kinks' or dispersion anoma- 
lies in the 50-70 meV energy range, which arise from the 
coupling of the electronic system with phonons|2| and/or 
magnetic modes[3]. Interestingly, the 'spaghetti' of vari- 
ous hybridized and unhybridized Cu and O bands start- 
ing around 1 eV in the cuprates seems to remain more or 
less unrenormalized. Q Recalling that opticaljH and elec- 
tron energy loss spectroscopy (EELS) experiments [l, 0] 
have shown the presence of dispersive plasmons in the 
cuprates lying at ~1 eV, we explore the effects of such a 
low energy plasmon in this Letter. The plasmon is found 
to induce renormalizations in the low energy regime in 
both the metallic and the insulating cuprates, which are 
in substantial accord with the corresponding experimen- 
tal results. We emphasize that plasmons in 3D materials 
typically lie at ~10 eV and therefore have little influence 
on the electronic states in the 0-1 eV range. 

We approach the problem within the framework of 
a one-band tight-binding model and consider the lim- 
iting cases of a metal with the example of overdoped 
Bi 2 Sr 2 CaCu20 8 (Bi2212)(lQ and the insulator with the 
example of half-filled Sr 2 Cu0 2 Cl 2 (SCOC) IS . The di- 
electric functions computed via the random phase ap- 
proximation (RPA) yield loss functions and the asso- 
ciated plasmon dispersions in both Bi2212 and SCOC, 
which are in reasonable accord with the corresponding 
EELS data. The generic effects of the electron-plasmon 
interaction on the electronic structure and how the low 
energy plasmon in the metal and the insulator induces 
characteristic dispersion anomalies and renormalizations 
of the spectrum are then delineated via a first order cal- 
culation of the electronic self energy. 

We discuss first the metallic case of overdoped Bi2212. 
Here the magnetic gap has collapsed fully and we take 
the bare electronic dispersion to be given by a one-band 
tight-binding model fitted to the LDA-based band struc- 
ture of Bi2212 in the vicinity of the Fermi energy (Ep) 



as follows^: 

e k = -2t[c x (a) + Cy(a)} - At' c x {a)c y {a) - 2t"[c x (2a) 

+c y (2a)] - At"'[c x (2a)c y (a) + c y {2a)c x {a)] , (1) 

where c a (na) = cos(nk a a) for a = x or y, n is an in- 
teger, and a is the in-plane lattice constant. Hopping 
parameters appropriate for Bi2212 (neglecting bilayer 
splitting) Q are: £=360 meV, t' = -100 mcV, t"=35 
meV, t'" = lOmeV. We model the dynamic dielectric 
function, e(q,(j), in the form (til. Hl| 
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where eo = 6 and the dispersion and damping of the 
plasmon in Bi2212 obtained from EELS data are given 



l.leV[l + AS xy ], 7 = 0.7eV[l + 16S ; 
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by[6j: uj p 

S xy = sin 2 (q x a/2) + sin 2 (q y a/2). Notably, due to the 
large broadening, the experimental plasmon peak lies at 
a lower energy than the bare value u> p . Using the bare 
dispersion of Eq. 1 and the dielectric function of Eq. 2, 
it is straightforward [12] to obtain the first order correc- 
tion to the electronic self energy S due to the screened 
Coulomb interaction [l3j. The formulas can also be de- 
rived from those for the insulator (Eqs. 5-11 below) by 
letting the magnetic gap A m — > 0. 

Figure 1 displays the spectral density in the overdoped 
metallic case where effects of the electron-plasmon inter- 
action are included. We consider the low energy region 
along r — > (it, 0) shown in (b) first. A region of high 
spectral intensity depicted by the thin white trace in (b) 
can be seen clearly (note logarithmic scale), indicating 
that the LDA band continues to persist in the presence 
of the plasmon, albeit with some renormalization. An 
average factor of Z = 0.5 is seen to more or less describe 
the renormalization of the LDA band as seen by compar- 
ing the dashed blue line with the whitish trace of high 
intensity. The Z = 0.5 value so determined is in substan- 
tial accord with the corresponding experimental value of 
Z = 0.28,[l| indicating that the plasmon constitutes a 
substantial part of the bosonic dressing of the electron. 
Moreover, the spectral weight spreads out from the renor- 
malized band towards the bare band (red line) , giving the 
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FIG. 1: (Color online) (a) Dressed spectral density function 
in overdoped (metallic) Bi2212. Bare dispersion (red line) is 
compared with dispersion renormalized with Z = 0.5 (blue 
dashed line). Dotted lines indicate plasmon dispersion, (b) 
Blowup of region enclosed by white box in (a). Logarithmic 
color scale is used to reveal weaker features. 



resulting yellowish trace the visual appearance of a 'wa- 
terfall' terminating around the T point. These results 
bear some resemblance to 'waterfall'-like anomalies re- 
ported recently in the ARPES spectra from the cuprates 
over the 0.4-0.8 eV energies However, the waterfall- 
like features induced by the plasmon and those reported 
in ARPES differ significantly insofar as their location in 
energy and momentum is concerned, and therefore, on 
the face of it, the plasmon would not seem to provide an 
explanation of the recent ARPES results. We return to 
comment further on this point below. 

We now turn to Fig. 1(a) which exposes effects of the 
plasmon over a wider energy scale. Here we see that 
the plasmon induces 'plasmaron' features 1J| - a new 



'shadow' band lying 1.5 cV below the bare dispersion via 
the plasmon emission process along the filled portion of 
the LDA bands along the T to (tt, 0) and the T to (tt, tt) 
lines. A similar shadow band associated with plasmon 
absorption appears 1.5 eV above the unfilled part of the 
LDA band along the (tt, 0) to (tt, tt) and T to (tt, tt) lines. 
At low energies, the plasmon dresses the electrons, mak- 
ing them heavier and renormalizing their dispersion. In 
contrast, at high energies, the plasmon becomes progres- 
sively less effective in dressing the electrons which move 
with ever increasing speeds, even though the electronic 
states are still broadened by plasmon emission effects. 

We now demonstrate that the experimentally observed 
plasmon dispersion in Bi2212 is in fact described reason- 
ably well by the band structure of Eq. 1, consistent with 
the results of Ref. [lj|. For this purpose, we directly 
evaluate the dielectric function from the bare charge sus- 
ceptibility xo an d the Coulomb interaction Vq via 



e(q,w) = e [l + V q xo(q,w)]. 



(3) 
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FIG. 2: (color online) Different calculations of the Coulomb 
potential (see text) in the form of qVq are compared. Full 3D 
layered lattice (solid line); A single 2D layer (dashed); and 
2D layer without the on-site U term (dotted). 



malized by a factor of Z=0.28 to account for the renor- 
malization observed in ARPES experiments on Bi2212. 
The treatment of the long-range Coulomb interaction 
V(r) for the correlated electronic system in a layered 
structure requires some care and we model V(r) as an 
on-site Hubbard U =2 cV term plus a 1/r contribu- 
tion screened by the background dielectric constant eo 
at other lattice sites. V q is then obtained via a lattice 
Fourier transform by summing the contributions of all 
in-plane Cu terms for r < R = 4a (a being the in-plane 
lattice constant) and approximating the remaining in- 
plane as well as interplane contributions for r > R by a 
continuum in order to recover the correct q — > limit of 



Xo is calculated straightforwardly within the RPA ap- 
proximation where the band structure of Eq. 1 is renor- 



V q . Fig. [2] shows qVq along high symmetry lines in the 
Brillouin zone (BZ). is seen to display a crossover from 
3D behavior ~ 1/q 2 to 2D behavior ~ 1/q for qc > 1, 
where c is the distance between successive Cu02 planes. 

The use of V q and \o in Eq. 3 yields the dielectric func- 
tion based directly on the band structure of Eq. 1 . The 
plasmon dispersion is found from the peak of the EELS 
spectrum, which is proportional to the loss function Im 
[— 1/e], plotted in Fig. [3] The theoretically predicted 
plasmon energy given by the maximum of the loss func- 
tion (white trace in the figure) is compared to the avail- 
able experimental results in Bi2212 (red dots). Given 
that the calculation is not self-consistent, the agreement 
is seen to be reasonable. Notably, our calculations as- 
sume purely in-plane excitations (q z = 0) . [17|| 

We consider next the more challenging case of insulat- 
ing SCOC, where a magnetic gap ispresent and it has 
not been clear how the EELS spectra[7( are related to the 
plasmon excitation. As in the metallic case, the starting 
point is the dispersion of Eq. 1, except that the hop- 
ping parameters are chosen to fit the LDA-based bands in 
SCOC: i=420 mcV, t' = -57 meV, £"=52 mcV, t'" = 26 
meV. The correlated band structure is modeled by intro- 
ducing an on-site Hubbard term U—2 eV and trea ting 
the system as a uniform half-filled antiferromagnet [18j 
with saturated (staggered) moment toq = 0.5. The band 
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FIG. 3: (Color online) Loss function Im [— 1/e] as a function 
of uj and q in Bi2212 plotted on a logarithmic color scale. Red 
dots give the experimental plasmon peak positions 



structure of Eq. 1 then splits into upper and lower mag- 



Um 



Q- 



The dielectric re- 



netic bands with gap A„ 
sponse given by Eq. 03 in the metallic case is modified 
strongly at low frequencies and our analysis indicates 
that the dielectric function in the insulator can be mod- 
eled as 



(4) 



where ujq ~ [/.[19j Taking oj =2 eV and the same 7 as 
in Eq. EH the EELS dataQ on SCOC can be fitted to 
Eq. |4]with Up = 1.3eF[l + 3S xy }. Again, the broadening 
shifts the peak in Im [e _1 ] from its expected position 
ujp = [ujp + ujq} 1 / 2 ; the experimental peak can be fit to a 
similar q-dependence, but with ujp = 2AeV[l + 0.335a;,;]. 
Note that ui p is not very different in metallic Bi2212 and 
insulating SCOC, so that most of the shift in the position 
u) p of the plasmon peak in the loss function is due to the 
gap w - 

The spectral density for the dressed Green function in 
the insulator is shown in Fig. 4. The results are similar 
to those of Fig. 1 for Bi2212 in that the dressed lower 
magnetic band in Fig. 4 displays the waterfall effect and 
that it is renormalized in comparison to the bare band 
(compare red line with dashed blue line). The renormal- 
ization of the low energy spectrum by a factor of Z ~ 0.3 
in SCOC due to the plasmon in Fig. 4 is again compara- 
ble to the corresponding experimental value of Z = 0.64. 

Finally, we consider briefly the calculation of the loss 
function and the self-energy in the insulating system. 
With two bands involved in this case, the formalism is 
less familiar than the one-band metallic case, and there- 
fore, an overview of the relvant results is appropriate. 
[The metallic case can be recovered by taking the limit 
A m — > 0.] The RPA expression for the charge suscepti- 
bility xo i n the insulator is 
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FIG. 4: (Color online) Dressed spectral density function in 
SCOC at half-filling. Bare dispersion (red line) is compared 
with dispersion renormalized with Z = 0.3 (blue dashed line). 
Dotted lines indicate plasmon dispersion. Logarithmic color 
scale is the same as in Fig. 1. 




(*,0) 



FIG. 5: (Color online) Loss function Im [—1/e] as a function 
of uj and q in SCOC at half-filling. Solid line gives the exper- 
imental plasmon dispersion (uj p ) for SCOC [7|. Logarithmic 
scale as in Fig. 3. 



where the prime means summation over the magnetic 
Brillouin zone only, 



2 = [1 ~ f(^i,k)]f{ej,k+q) 

9ij 



£i,k - £j,k+q + w + i6' 
f is the Fermi function, 



£l,2fe 



-,(+) 



with the 1(2) going with the -(+) sign, 
£fc+<?)/2, 



1 r 
2 



1 ± 



e (-) e (-) 



E k E k 



+q 



(6) 

(7) 
(e fc ± 

(8) 



and Ek = [e k ^ 2 + A^J 1 / 2 . The preceding equations have 
been applied previously to the Hubbard model [13, 21 1, 
but we have extended their usage to include the long- 
range Coulomb interaction. 
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In the presence of a gap, the electronic self energy £ 
becomes a tensor with indices 1(2) corresponding to the 



lower (upper)magnetic band[12l. [2 



S u (p,w,o-) 



BZ 



S12 — £21 — 



BZ 



(2tt)* 



(9) 



(10) 



where 



= [Q(oj - e liP+q ) - Q(-ei tP+q )] 




^ q 2 Xo(q^ ~ £i,p+q) _ 
e(q, oj - ei, p+g )/eo ' 

e (-) 

P+9 




E 



p+q 



(11) 



and 



and £22 has the same form as En, but with u n 
interchanged. 

The loss function in SCOC computed from Eqs. 5 and 
3 is shown in Fig. 5. The experimentally observed plas- 
mon dispersion (solid line) is seen to be in good accord 
with the theoretical prediction. A comparison of Figs. 
3 and 5 reveals how the gap in the underlying spectrum 
of the insulator shifts the plasmon energy in SCOC to 
higher values with little weight at energies below w 2 eV 
in contrast to the metallic case of Bi2212. 

The present results have significant implications for the 
low energy physics of the cuprates and provide a route 
for connecting this physics to the spectrum at higher en- 
ergy scales as one encounters the venerable 'spaghetti' of 
Cu and O related bands spread typically over ~l-8 eV 
binding energies. We have shown that the plasmon in- 
duces band renormalizations at low energies in the metal- 
lic case of Bi2212 as well as the insulating case of SCOC, 
and that these renormalization effects are in substantial 
accord with the corresponding experimental values. Our 
results clearly establish the importance of the plasmon 
in the bosonic dressings of the low energy quasiparticles 
in the cuprates. 

There have been several recent calculations of 
electronically-induced kinks, but these are gene rally very 
wcak|23] or confined to a restricted doping range 24]. The 
plasmon by contrast is robust and produces strong kink- 
like features which extend over several eVs, due to the 
large broadening of the plasmon 25j • While these char- 
acteristic features of the electron-plasmon interaction are 
reminiscent of those reported in very recent ARPES ex- 
periments on several families of cuprates [H, 0, Eo| . they 
fall at too high an energy scale. Possible candidates 
for the intermediate waterfall scale include the acoustic 
branch of the plasmons and high-energy magnetic exci- 
tations. The extent to which low energy plasmons are 
involved in an electronic mechanism of superconductiv- 
ity remains unclear, although such suggestions have been 
made in the literature 12611. 



In conclusion, we have shown that the low-energy plas- 
mon, which is a universal feature of the cuprates, pro- 
duces renormalizations and dispersion anomalies in the 
electronic structure. Our study indicates that the plas- 
mon is a very significant player for understanding the 
physics of the cuprates. 
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